Novel signals and polygenic score for height are associated with pubertal growth traits in Southwestern American Indians.

Most genetic variants associated with adult height have been identified through large genome-wide association studies (GWASs) in European-ancestry cohorts. However, it is unclear how these variants influence linear growth during adolescence. This study uses anthropometric and genotypic data from a longitudinal study conducted in an American Indian community in Arizona between 1965-2007. Growth parameters (i.e. height, velocity, and timing of growth spurt) were derived from the Preece-Baines growth model, a parametric growth curve fitted to longitudinal height data, in 787 participants with height measurements spanning the whole period of growth. Heritability estimates suggested that genetic factors could explain 25% to 71% of the variance of pubertal growth traits. We performed a GWAS of growth parameters, testing their associations with 5 077 595 imputed or directly genotyped variants. Six variants associated with height at peak velocity (P < 5 × 10-8, adjusted for sex, birth year and principal components). Implicated genes include NUDT3, previously associated with adult height, and PACSIN1. Two novel variants associated with duration of growth spurt (P < 5 × 10-8) in LOC105375344, an uncharacterized gene with unknown function. We finally examined the association of growth parameters with a polygenic score for height derived from 9557 single nucleotide polymorphisms (SNPs) identified in the GIANT meta-analysis for which genotypic data were available for the American Indian study population. Height polygenic score was correlated with the magnitude and velocity of height growth that occurred before and at the peak of the adolescent growth spurt, indicating overlapping genetic architecture, with no influence on the timing of adolescent growth.

The utility of a type 2 diabetes polygenic score in addition to clinical variables for prediction of type 2 diabetes incidence in birth, youth and adult cohorts in an Indigenous study population.

AIMS/HYPOTHESIS: There is limited information on how polygenic scores (PSs), based on variants from genome-wide association studies (GWASs) of type 2 diabetes, add to clinical variables in predicting type 2 diabetes incidence, particularly in non-European-ancestry populations. METHODS: For participants in a longitudinal study in an Indigenous population from the Southwestern USA with high type 2 diabetes prevalence, we analysed ten constructions of PS using publicly available GWAS summary statistics. Type 2 diabetes incidence was examined in three cohorts of individuals without diabetes at baseline. The adult cohort, 2333 participants followed from age ≥20 years, had 640 type 2 diabetes cases. The youth cohort included 2229 participants followed from age 5-19 years (228 cases). The birth cohort included 2894 participants followed from birth (438 cases). We assessed contributions of PSs and clinical variables in predicting type 2 diabetes incidence. RESULTS: Of the ten PS constructions, a PS using 293 genome-wide significant variants from a large type 2 diabetes GWAS meta-analysis in European-ancestry populations performed best. In the adult cohort, the AUC of the receiver operating characteristic curve for clinical variables for prediction of incident type 2 diabetes was 0.728; with the PS, 0.735. The PS's HR was 1.27 per SD (p=1.6 × 10-8; 95% CI 1.17, 1.38). In youth, corresponding AUCs were 0.805 and 0.812, with HR 1.49 (p=4.3 × 10-8; 95% CI 1.29, 1.72). In the birth cohort, AUCs were 0.614 and 0.685, with HR 1.48 (p=2.8 × 10-16; 95% CI 1.35, 1.63). To further assess the potential impact of including PS for assessing individual risk, net reclassification improvement (NRI) was calculated: NRI for the PS was 0.270, 0.268 and 0.362 for adult, youth and birth cohorts, respectively. For comparison, NRI for HbA1c was 0.267 and 0.173 for adult and youth cohorts, respectively. In decision curve analyses across all cohorts, the net benefit of including the PS in addition to clinical variables was most pronounced at moderately stringent threshold probability values for instituting a preventive intervention. CONCLUSIONS/INTERPRETATION: This study demonstrates that a European-derived PS contributes significantly to prediction of type 2 diabetes incidence in addition to information provided by clinical variables in this Indigenous study population. Discriminatory power of the PS was similar to that of other commonly measured clinical variables (e.g. HbA1c). Including type 2 diabetes PS in addition to clinical variables may be clinically beneficial for identifying individuals at higher risk for the disease, especially at younger ages.

Identification and characterization of the long non-coding RNA NFIA-AS2 as a novel locus for body mass index in American Indians.

BACKGROUND: Genome-wide association studies have shown that body mass index (BMI), an estimate of obesity, is highly polygenic. Individual variants typically have small effect sizes, making it challenging to identify unique loci in under-represented ethnic groups which lack statistical power due to their small sample size. Yet obesity is a major health disparity and is particularly prevalent in southwestern American Indians. Here, we identify and characterize a new locus for BMI that was detected by analyzing moderate associations with BMI obtained in a population-based sample of southwestern American Indians together with the well-powered GIANT dataset. METHODS: Genotypes for 10.5 million variants were tested for association with BMI in 5870 American Indians and 2600 variants that showed an association P < 10-3 in the American Indian sample were combined in a meta-analysis with the BMI data reported in GIANT (N = 240,608). The newly identified gene, NFIA-AS2 was functionally characterized, and the impact of its lead associated variant rs1777538 was studied both in-silico and in-vitro. RESULTS: Rs1777538 (T/C; C allele frequency = 0.16 in American Indians and 0.04 in GIANT, meta-analysis P = 5.0 × 10-7) exhibited a large effect in American Indians (1 kg/m2 decrease in BMI per copy of C allele). NFIA-AS2 was found to be a nuclear localized long non-coding RNA expressed in tissues pertinent to human obesity. Analysis of this variant in human brown preadipocytes showed that NFIA-AS2 transcripts carrying the C allele had increased RNA degradation compared to the T allele transcripts (half-lives = 9 h, 13 h respectively). During brown adipogenesis, NFIA-AS2 featured a stage-specific regulation of nearby gene expression where rs1777538 demonstrated an allelic difference in regulation in the mature adipocytes (the strongest difference was observed for L1TD1, P = 0.007). CONCLUSION: Our findings support a role for NFIA-AS2 in regulating pathways that impact BMI.

Further evidence supporting a potential role for ADH1B in obesity.

Insulin is an essential hormone that regulates glucose homeostasis and metabolism. Insulin resistance (IR) arises when tissues fail to respond to insulin, and it leads to serious health problems including Type 2 Diabetes (T2D). Obesity is a major contributor to the development of IR and T2D. We previously showed that gene expression of alcohol dehydrogenase 1B (ADH1B) was inversely correlated with obesity and IR in subcutaneous adipose tissue of Mexican Americans. In the current study, a meta-analysis of the relationship between ADH1B expression and BMI in Mexican Americans, African Americans, Europeans, and Pima Indians verified that BMI was increased with decreased ADH1B expression. Using established human subcutaneous pre-adipocyte cell lines derived from lean (BMI < 30 kg m-2) or obese (BMI ≥ 30 kg m-2) donors, we found that ADH1B protein expression increased substantially during differentiation, and overexpression of ADH1B inhibited fatty acid binding protein expression. Mature adipocytes from lean donors expressed ADH1B at higher levels than obese donors. Insulin further induced ADH1B protein expression as well as enzyme activity. Knockdown of ADH1B expression decreased insulin-stimulated glucose uptake. Our findings suggest that ADH1B is involved in the proper development and metabolic activity of adipose tissues and this function is suppressed by obesity.

Assessment of the potential role of natural selection in type 2 diabetes and related traits across human continental ancestry groups: comparison of phenotypic with genotypic divergence.

AIMS/HYPOTHESIS: Prevalence of type 2 diabetes differs among human ancestry groups, and many hypotheses invoke differential natural selection to account for these differences. We sought to assess the potential role of differential natural selection across major continental ancestry groups for diabetes and related traits, by comparison of genetic and phenotypic differences. METHODS: This was a cross-sectional comparison among 734 individuals from an urban sample (none of whom was more closely related to another than third-degree relatives), including 83 African Americans, 523 American Indians and 128 European Americans. Participants were not recruited based on diabetes status or other traits. BMI was calculated, and diabetes was diagnosed by a 75 g oral glucose tolerance test. In those with normal glucose tolerance (n = 434), fasting insulin and 30 min post-load insulin, adjusted for 30 min glucose, were taken as measures of insulin resistance and secretion, respectively. Whole exome sequencing was performed, resulting in 97,388 common (minor allele frequency ≥ 5%) variants; the coancestry coefficient (FST) was calculated across all markers as a measure of genetic divergence among ancestry groups. The phenotypic divergence index (PST) was also calculated from the phenotypic differences and heritability (which was estimated from genetic relatedness calculated empirically across all markers in 761 American Indian participants prior to the exclusion of close relatives). Under evolutionary neutrality, the expectation is PST = FST, while for traits under differential selection PST is expected to be significantly greater than FST. A bootstrap procedure was used to test the hypothesis PST = FST. RESULTS: With adjustment for age and sex, prevalence of type 2 diabetes was 34.0% in American Indians, 12.4% in African Americans and 10.4% in European Americans (p = 2.9 × 10-10 for difference among groups). Mean BMI was 36.3, 33.4 and 33.0 kg/m2, respectively (p = 1.9 × 10-7). Mean fasting insulin was 63.8, 48.4 and 45.2 pmol/l (p = 9.2 × 10-5), while mean 30 min insulin was 559.8, 553.5 and 358.8 pmol/l, respectively (p = 5.7 × 10-8). FST across all markers was 0.130, while PST for liability to diabetes, adjusted for age and sex, was 0.149 (p = 0.35 for difference with FST). PST was 0.094 for BMI (p = 0.54), 0.095 for fasting insulin (p = 0.54) and 0.216 (p = 0.18) for 30 min insulin. For type 2 diabetes and BMI, the maximum divergence between populations was observed between American Indians and European Americans (PST-MAX = 0.22, p = 0.37, and PST-MAX = 0.14, p = 0.61), which suggests that a relatively modest 22% or 14% of the genetic variance, respectively, can potentially be explained by differential selection (assuming the absence of neutral drift). CONCLUSIONS/INTERPRETATION: These analyses suggest that while type 2 diabetes and related traits differ significantly among continental ancestry groups, the differences are consistent with neutral expectations based on heritability and genetic distances. While these analyses do not exclude a modest role for natural selection, they do not support the hypothesis that differential natural selection is necessary to explain the phenotypic differences among these ancestry groups. Graphical abstract.

Exome Sequencing Identifies A Nonsense Variant in DAO Associated With Reduced Energy Expenditure in American Indians.

BACKGROUND: Obesity and energy expenditure (EE) are heritable and genetic variants influencing EE may contribute to the development of obesity. We sought to identify genetic variants that affect EE in American Indians, an ethnic group with high prevalence of obesity. METHODS: Whole-exome sequencing was performed in 373 healthy Pima Indians informative for 24-hour EE during energy balance. Genetic association analyses of all high-quality exonic variants (≥5 carriers) was performed, and those predicted to be damaging were prioritized. RESULTS: Rs752074397 introduces a premature stop codon (Cys264Ter) in DAO and demonstrated the strongest association for 24-hour EE, where the Ter allele associated with substantially lower 24-hour EE (mean lower by 268 kcal/d) and sleeping EE (by 135 kcal/d). The Ter allele has a frequency = 0.5% in Pima Indians, whereas is extremely rare in most other ethnic groups (frequency < 0.01%). In vitro functional analysis showed reduced protein levels for the truncated form of DAO consistent with increased protein degradation. DAO encodes D-amino acid oxidase, which is involved in dopamine synthesis which might explain its role in modulating EE. CONCLUSION: Our results indicate that a nonsense mutation in DAO may influence EE in American Indians. Identification of variants that influence energy metabolism may lead to new pathways to treat human obesity. CLINICAL TRIAL REGISTRATION NUMBER: NCT00340132.

Characterization of Exome Variants and Their Metabolic Impact in 6,716 American Indians from the Southwest US.

Applying exome sequencing to populations with unique genetic architecture has the potential to reveal novel genes and variants associated with traits and diseases. We sequenced and analyzed the exomes of 6,716 individuals from a Southwestern American Indian (SWAI) population with well-characterized metabolic traits. We found that the SWAI population has distinct allelic architecture compared to populations of European and East Asian ancestry, and there were many predicted loss-of-function (pLOF) and nonsynonymous variants that were highly enriched or private in the SWAI population. We used pLOF and nonsynonymous variants in the SWAI population to evaluate gene-burden associations of candidate genes from European genome-wide association studies (GWASs) for type 2 diabetes, body mass index, and four major plasma lipids. We found 19 significant gene-burden associations for 11 genes, providing additional evidence for prioritizing candidate effector genes of GWAS signals. Interestingly, these associations were mainly driven by pLOF and nonsynonymous variants that are unique or highly enriched in the SWAI population. Particularly, we found four pLOF or nonsynonymous variants in APOB, APOE, PCSK9, and TM6SF2 that are private or enriched in the SWAI population and associated with low-density lipoprotein (LDL) cholesterol levels. Their large estimated effects on LDL cholesterol levels suggest strong impacts on protein function and potential clinical implications of these variants in cardiovascular health. In summary, our study illustrates the utility and potential of exome sequencing in genetically unique populations, such as the SWAI population, to prioritize candidate effector genes within GWAS loci and to find additional variants in known disease genes with potential clinical impact.

Low Serum Insulinlike Growth Factor II Levels Correlate with High BMI in American Indian Adults.

OBJECTIVE: Insulinlike growth factor II (IGF-II) regulates metabolism and growth. In humans, both positive and negative relationships have been reported between serum IGF-II levels and obesity. This study assessed the relationship between serum IGF-II levels and BMI and determined whether IGF-II levels predict weight gain. METHODS: Serum samples were available from 911 American Indians with a recorded BMI. IGF-II was measured using enzyme-linked immunosorbent assay. RESULTS: Serum IGF-II levels were negatively correlated with BMI (r = -0.17, P = 4.4 × 10-7 , adjusted for age, sex, and storage time). The strongest correlation was in participants aged ≥ 30 years (r = -0.28, P = 3.4 × 10-8 , N = 349), a modest correlation was in participants aged 20 to 29 years (r = -0.15, P = 7.6 × 10-3 , N = 322), and participants aged 15 to 19 years had no correlation (r = 0.05, P = 0.48, N = 240). IGF-II levels did not predict weight gain. However, among individuals who had genotypes for 64 established obesity variants (age ≥ 20 years, N = 671), a genetic risk score for high BMI was associated with lower IGF-II (ß = -0.08 SD of IGF-II per SD of the genetic risk score, P = 0.025). CONCLUSIONS: There is a negative relationship between IGF-II levels and BMI, in which the correlation is stronger at older ages. The association between genetic risk for BMI and IGF-II levels suggests that this correlation may be due to an effect of obesity on IGF-II.

Association of CREBRF variants with obesity and diabetes in Pacific Islanders from Guam and Saipan.

AIMS/HYPOTHESIS: Variants in CREBRF (rs12513649 and rs373863828) have been strongly associated with increased BMI and decreased risk of type 2 diabetes in Polynesian populations; the A allele at rs373863828 is common in Polynesians but rare in most other global populations. The aim of the present study was to assess the association of CREBRF variants with obesity and diabetes in Pacific Islander (largely Marianas and Micronesian) populations from Guam and Saipan. METHODS: CREBRF rs12513649 and rs373863828 were genotyped in 2022 participants in a community-based cross-sectional study designed to identify determinants of diabetes and end-stage renal disease (ESRD). Associations were analysed with adjustment for age, sex, ESRD and the first four genetic principal components from a genome-wide association study (to account for population stratification); a genomic control procedure was used to account for residual stratification. RESULTS: The G allele at rs12513649 had an overall frequency of 7.7%, which varied from 2.2% to 20.7% across different Marianas and Micronesian populations; overall frequency of the A allele at rs373863828 was 4.2% (range: 1.1-5.4%). The G allele at rs12513649 was associated with higher BMI (ß = 1.55 kg/m2 per copy; p = 0.0026) as was the A allele at rs373863828 (ß = 1.48 kg/m2, p = 0.033). The same alleles were associated with lower risk of diabetes (OR per copy: 0.63 [p = 0.0063] and 0.49 [p = 0.0022], respectively). Meta-analyses combining the current results with previous results in Polynesians showed a strong association between the A allele at rs373863828 and BMI (ß = 1.38 kg/m2; p = 2.5 × 10-29) and diabetes (OR 0.65, p = 1.5 × 10-13). CONCLUSIONS/INTERPRETATION: These results confirm the associations of CREBRF variants with higher BMI and lower risk of diabetes and, importantly, they suggest that these variants contribute to the risk of obesity and diabetes in Oceanic populations.

Assessing the Role of 98 Established Loci for BMI in American Indians.

OBJECTIVE: Meta-analyses of genome-wide association studies in Europeans have identified > 98 loci for BMI. Transferability of these established associations in Pima Indians was analyzed. METHODS: Among 98 lead single nucleotide polymorphisms (SNPs), 82 had minor allele frequency ≥ 0.01 in Pima Indians and were analyzed for association with the maximum BMI in adulthood (n = 3,491) and BMI z score in childhood (n = 1,958). Common tag SNPs across 98 loci were also analyzed for additional signals. RESULTS: Among the lead SNPs, 13 (TMEM18, TCF7L2, MRPS33P4, PRKD1, ZFP64, FTO, TAL1, CALCR, GNPDA2, CREB1, LMX1B, ADCY9, NLRC3) were associated with BMI (P  ≤ 0.05) in Pima adults. A multi-allelic genetic risk score (GRS), which summed the risk alleles for 82 lead SNPs, showed a significant trend for a positive relationship between GRS and BMI in adulthood (beta = 0.48% per risk allele; P = 1.6 × 10-9 ) and BMI z score in childhood (beta = 0.024 SD; P = 1.7 × 10-7 ). GRS was significantly associated with BMI across all age groups ≥ 5 years, except for those ≥ 50 years. The strongest association was seen in adolescence (age 14-16 years; P = 1.84 × 10-9 ). CONCLUSIONS: In aggregate, European-derived lead SNPs had a notable effect on BMI in Pima Indians. Polygenic obesity in this population manifests strongly in childhood and adolescence and persists throughout much of adult life.

Analysis of type 2 diabetes and obesity genetic variants in Mexican Pima Indians: Marked allelic differentiation among Amerindians at HLA.

Prevalence of diabetes and obesity in Mexican Pima Indians is low, while prevalence in US Pima Indians is high. Although lifestyle likely accounts for much of the difference, the role of genetic factors is not well explored. To examine this, we genotyped 359 single nucleotide polymorphisms, including established type 2 diabetes and obesity variants from genome-wide association studies (GWAS) and 96 random markers, in 342 Mexican Pimas. A multimarker risk score of obesity variants was associated with body mass index (BMI; ß = 0.81 kg/m2 per SD, P = 0.0066). The mean value of the score was lower in Mexican Pimas than in US Pimas (P = 4.3 × 10-11 ), and differences in allele frequencies at established loci could account for approximately 7% of the population difference in BMI; however, the difference in risk scores was consistent with evolutionary neutrality given genetic distance. To identify loci potentially under recent natural selection, allele frequencies at 283 variants were compared between US and Mexican Pimas, accounting for genetic distance. The largest differences were seen at HLA markers (e.g., rs9271720, difference = 0.75, P = 8.7 × 10-9 ); genetic distances at HLA were greater than at random markers (P = 1.6 × 10-46 ). Analyses of GWAS data in 937 US Pimas also showed sharing of alleles identical by descent at HLA that exceeds its genomic expectation (P = 7.0 × 10-10 ). These results suggest that, in addition to the widely recognized balancing selection at HLA, recent directional selection may also occur, resulting in marked allelic differentiation between closely related populations.

Identification and functional analysis of a novel G310D variant in the insulin-like growth factor 1 receptor (IGF1R) gene associated with type 2 diabetes in American Indians.

AIMS: Insulin-like growth factor 1 receptor (IGF1R) is involved in cell growth and glucose homeostasis. In the current study, the IGF1R locus was analysed as a candidate gene for type 2 diabetes (T2D) in American Indians. MATERIALS AND METHODS: Whole genome sequence data from 335 American Indians identified 3 novel missense variants in IGF1R. The associations of IGF1R variants with T2D, age of T2D onset and birth weight were analysed in a population-based sample of 7701 American Indians. RESULTS: A novel glycine-to-aspartic acid substitution (G310D) in IGF1R was identified, which associated with T2D in a sex-specific manner (Psex interaction = 0.02). In women, the aspartic acid (D) allele (frequency = 0.034) was associated with increased risk for T2D (n = 4292, P = 2.0 × 10-5 adjusted for age, birth year, and the first 5 genetic principal components; odds ratio [OR] = 2.23 [1.54-3.23] per risk allele) and an earlier age of T2D onset (n = 4292, P = 2 × 10-4 , hazard rate ratio = 1.45 [1.20-1.75], Psex interaction = 0.05). Female carriers of the D-allele also had lower birth weight (n = 1313, ß = -163 g, P = .006, Psex interaction = 0.008). Among 85 siblings discordant for G310D, carriers of the D-allele had shorter stature as compared with carriers of the G-allele (ß = -1.6 cm, P = .001, within family model). The G310D variant was functionally studied in vitro, where the D-allele had a 22% increase (P = .0005) in FOXO1-induced transcriptional activity, due to decreased activation of the PI3K/AKT pathway mediated through reduced IGF1R activity. CONCLUSION: A unique G310D variant in IGF1R, which occurs in 6% American Indians, may impair IGF1R signalling pathways, thereby increasing the risk of T2D.

Serum lipids and mortality in an American Indian population: A longitudinal study.

BACKGROUND: In Caucasians, lower triglycerides (TG), total or LDL cholesterol and high HDL cholesterol are generally associated with lower mortality. However, low cholesterol is associated with higher mortality in some Asian populations. This study examines the relationship between serum lipids and mortality in American Indians. METHODS: 2125 American Indians aged ≥40years were examined biennially between 1993 and 2007. Vital status was determined through 2011. Mortality rates, adjusted for age, sex and diabetes, were calculated using Poisson regression. RESULTS: The median baseline age was 46years and 61% were women. Over a median follow-up of 10.1years, 522 deaths occurred. Relationships between baseline lipids, except for HDL cholesterol, and all-cause mortality were negative and linear in persons without diabetes and U-shaped in persons with diabetes. For HDL cholesterol, the relationship was U-shaped in the total cohort. Cardiovascular mortality was positively associated with total, LDL and non-HDL cholesterol whereas lower lipid concentrations were adversely associated with mortality from liver disease or external causes, except for HDL cholesterol, where associations were positive. CONCLUSION: The common belief that low cholesterol and TG are beneficial for health is not universally observed; evidence suggests increased mortality at both ends of the cholesterol and TG distributions.

Identity-by-Descent Mapping Identifies Major Locus for Serum Triglycerides in Amerindians Largely Explained by an APOC3 Founder Mutation.

BACKGROUND: Identity-by-descent mapping using empirical estimates of identity-by-descent allele sharing may be useful for studies of complex traits in founder populations, where hidden relationships may augment the inherent genetic information that can be used for localization. METHODS AND RESULTS: Through identity-by-descent mapping, using ≈400 000 single-nucleotide polymorphisms (SNPs), of serum lipid profiles, we identified a major linkage signal for triglycerides in 1007 Pima Indians (LOD=9.23; P=3.5×10-11 on chromosome 11q). In subsequent fine-mapping and replication association studies in ≈7500 Amerindians, we determined that this signal reflects effects of a loss-of-function Ala43Thr substitution in APOC3 (rs147210663) and 3 established functional SNPs in APOA5. The association with rs147210663 was particularly strong; each copy of the Thr allele conferred 42% lower triglycerides (ß=-0.92±0.059 SD unit; P=9.6×10-55 in 4668 Pimas and 2793 Southwest Amerindians combined). The Thr allele is extremely rare in most global populations but has a frequency of 2.5% in Pimas. We further demonstrated that 3 APOA5 SNPs with established functional impact could explain the association with the most well-replicated SNP (rs964184) for triglycerides identified by genome-wide association studies. Collectively, these 4 SNPs account for 6.9% of variation in triglycerides in Pimas (and 4.1% in Southwest Amerindians), and their inclusion in the original linkage model reduced the linkage signal to virtually null. CONCLUSIONS: APOC3/APOA5 constitutes a major locus for serum triglycerides in Amerindians, especially the Pimas, and these results provide an empirical example for the concept that population-based linkage analysis is a useful strategy to identify complex trait variants.

A Genome-Wide Association Study Using a Custom Genotyping Array Identifies Variants in GPR158 Associated With Reduced Energy Expenditure in American Indians.

Pima Indians living in Arizona have a high prevalence of obesity, and we have previously shown that a relatively lower energy expenditure (EE) predicts weight and fat mass gain in this population. EE is a familial trait (heritability = 0.52); therefore, in the current study, we aimed to identify genetic variants that affect EE and thereby influence BMI and body fatness in Pima Indians. Genotypic data from 491,265 variants were analyzed for association with resting metabolic rate (RMR) and 24-h EE assessed in a whole-room calorimeter in 507 and 419 Pima Indians, respectively. Variants associated with both measures of EE were analyzed for association with maximum BMI and percent body fat (PFAT) in 5,870 and 912 Pima Indians, respectively. rs11014566 nominally associated with both measures of EE and both measures of adiposity in Pima Indians, where the G allele (frequency: Pima Indians = 0.60, Europeans <0.01) associated with lower 24-h EE (ß = -33 kcal/day per copy), lower RMR (ß = -31 kcal/day), higher BMI (ß = +0.6 kg/m2), and higher PFAT (ß = +0.9%). However, the association of rs11014566 with BMI did not directionally replicate when assessed in other ethnic groups. rs11014566 tags rs144895904, which affected promoter function in an in vitro luciferase assay. These variants map to GPR158, which is highly expressed in the brain and interacts with two other genes (RGS7 and CACNA1B) known to affect obesity in knockout mice. Our results suggest that common ethnic-specific variation in GPR158 may influence EE; however, its role in weight gain remains controversial, as it either had no association with BMI or associated with BMI but in the opposite direction in other ethnic groups.

Fine-mapping, novel loci identification, and SNP association transferability in a genome-wide association study of QRS duration in African Americans.

The electrocardiographic QRS duration, a measure of ventricular depolarization and conduction, is associated with cardiovascular mortality. While single nucleotide polymorphisms (SNPs) associated with QRS duration have been identified at 22 loci in populations of European descent, the genetic architecture of QRS duration in non-European populations is largely unknown. We therefore performed a genome-wide association study (GWAS) meta-analysis of QRS duration in 13,031 African Americans from ten cohorts and a transethnic GWAS meta-analysis with additional results from populations of European descent. In the African American GWAS, a single genome-wide significant SNP association was identified (rs3922844, P = 4 × 10-14) in intron 16 of SCN5A, a voltage-gated cardiac sodium channel gene. The QRS-prolonging rs3922844 C allele was also associated with decreased SCN5A RNA expression in human atrial tissue (P = 1.1 × 10-4). High density genotyping revealed that the SCN5A association region in African Americans was confined to intron 16. Transethnic GWAS meta-analysis identified novel SNP associations on chromosome 18 in MYL12A (rs1662342, P = 4.9 × 10-8) and chromosome 1 near CD1E and SPTA1 (rs7547997, P = 7.9 × 10-9). The 22 QRS loci previously identified in populations of European descent were enriched for significant SNP associations with QRS duration in African Americans (P = 9.9 × 10-7), and index SNP associations in or near SCN5A, SCN10A, CDKN1A, NFIA, HAND1, TBX5 and SETBP1 replicated in African Americans. In summary, rs3922844 was associated with QRS duration and SCN5A expression, two novel QRS loci were identified using transethnic meta-analysis, and a significant proportion of QRS-SNP associations discovered in populations of European descent were transferable to African Americans when adequate power was achieved.

The Association of Serum Leptin with Mortality in Older Adults.

OBJECTIVE: Elevated levels of serum leptin are associated with increased adiposity and production of pro-inflammatory cytokines. Both cytokines and body adiposity have been shown to predict cardiovascular events and mortality. The primary objective of the present study is to explore the associations between serum leptin and all-cause mortality and mortality from cardiovascular disease (CVD) over a span of 10 years, controlling for body adiposity and proinflammatory cytokines. METHODS: The Health, Aging and Body Composition (Health ABC) study is a prospective cohort of 3,075 older adults aged 70 to 79 years. This analysis includes 2,919 men and women with complete serum leptin and vital status data. Data on all-cause mortality and incident cardiovascular events (including Coronary Heart Disease and Congestive Heart Failure) were collected over 10 years of follow-up (mean 8.4 years). RESULTS: Women with leptin in quartile 2 and 3 were at lower risk of all-cause mortality, and those with leptin in quartile 2 were at lower risk of mortality from CVD as compared to women with lowest leptin values when adjusted for age, race, site, years of education, alcohol use, smoking, and physical activity. When these associations were additionally adjusted for body fat, C-reactive protein and pro-inflammatory cytokines, women with leptin values in quartile 3 were at lower risk of all-cause mortality and women with leptin in quartile 2 and 3 were at lower risk of mortality from CVD than women with lowest leptin values. These associations were not significant among men after adjusting for body fat and cytokines. CONCLUSIONS: The present study suggests that moderately elevated concentrations of serum leptin are independently associated with lower risk of all-cause mortality and CVD-related mortality among older women. Among men, serum leptin is not associated with reduced risk of all-cause and CVD mortality after controlling for body fat and cytokines.

Chronotype and seasonality: morningness is associated with lower seasonal mood and behavior changes in the Old Order Amish.

BACKGROUND: Several studies documented that lower scores on the Morningness-Eveningness Questionnaire (MEQ) are associated with a higher global seasonality of mood (GSS). As for the Modern Man artificial lighting predominantly extends evening activity and exposure to light, and as evening bright light phase is known to delay circadian rhythms, this chronic exposure could potentially lead to both lower Morningness as well as higher GSS. The aim of the study was to investigate if the MEQ-GSS relationship holds in the Old Order Amish of Lancaster County, PA, a population that does not use network electrical light. METHODS: 489 Old Order Amish adults (47.6% women), with average (SD) age of 49.7 (14.2) years, completed both the Seasonal Pattern Assessment Questionnaire (SPAQ) for the assessment of GSS, and MEQ. Associations between GSS scores and MEQ scores were analyzed using linear models, accounting for age, gender and relatedness by including the relationship matrix in the model as a random effect. RESULTS: GSS was inversely associated with MEQ scores (p=0.006, adjusted). LIMITATIONS: include a potential recall bias associated with self-report questionnaires and no actual light exposure measurements. CONCLUSION: We confirmed the previously reported inverse association between MEQ scores and lower seasonality of mood, for the first time in a population that does not use home network electrical lighting. This result suggests that the association is not a byproduct of exposure to network electric light, and calls for additional research to investigate mechanisms by which Morningness is negatively associated with seasonality.

Genetic association study of adiposity and melanocortin-4 receptor (MC4R) common variants: replication and functional characterization of non-coding regions.

Common genetic variants 3' of MC4R within two large linkage disequilibrium (LD) blocks spanning 288 kb have been associated with common and rare forms of obesity. This large association region has not been refined and the relevant DNA segments within the association region have not been identified. In this study, we investigated whether common variants in the MC4R gene region were associated with adiposity-related traits in a biracial population-based study. Single nucleotide polymorphisms (SNPs) in the MC4R region were genotyped with a custom array and a genome-wide array and associations between SNPs and five adiposity-related traits were determined using race-stratified linear regression. Previously reported associations between lower BMI and the minor alleles of rs2229616/Val103Ile and rs52820871/Ile251Leu were replicated in white female participants. Among white participants, rs11152221 in a proximal 3' LD block (closer to MC4R) was significantly associated with multiple adiposity traits, but SNPs in a distal 3' LD block (farther from MC4R) were not. In a case-control study of severe obesity, rs11152221 was significantly associated. The association results directed our follow-up studies to the proximal LD block downstream of MC4R. By considering nucleotide conservation, the significance of association, and proximity to the MC4R gene, we identified a candidate MC4R regulatory region. This candidate region was sequenced in 20 individuals from a study of severe obesity in an attempt to identify additional variants, and the candidate region was tested for enhancer activity using in vivo enhancer assays in zebrafish and mice. Novel variants were not identified by sequencing and the candidate region did not drive reporter gene expression in zebrafish or mice. The identification of a putative insulator in this region could help to explain the challenges faced in this study and others to link SNPs associated with adiposity to altered MC4R expression.

Who are the Okinawans? Ancestry, genome diversity, and implications for the genetic study of human longevity from a geographically isolated population.

Isolated populations have advantages for genetic studies of longevity from decreased haplotype diversity and long-range linkage disequilibrium. This permits smaller sample sizes without loss of power, among other utilities. Little is known about the genome of the Okinawans, a potential population isolate, recognized for longevity. Therefore, we assessed genetic diversity, structure, and admixture in Okinawans, and compared this with Caucasians, Chinese, Japanese, and Africans from HapMap II, genotyped on the same Affymetrix GeneChip Human Mapping 500K array. Principal component analysis, haplotype coverage, and linkage disequilibrium decay revealed a distinct Okinawan genome-more homogeneity, less haplotype diversity, and longer range linkage disequilibrium. Population structure and admixture analyses utilizing 52 global reference populations from the Human Genome Diversity Cell Line Panel demonstrated that Okinawans clustered almost exclusively with East Asians. Sibling relative risk (λs) analysis revealed that siblings of Okinawan centenarians have 3.11 times (females) and 3.77 times (males) more likelihood of centenarianism. These findings suggest that Okinawans are genetically distinct and share several characteristics of a population isolate, which are prone to develop extreme phenotypes (eg, longevity) from genetic drift, natural selection, and population bottlenecks. These data support further exploration of genetic influence on longevity in the Okinawans.